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Abstract 

We present a picture of phase transitions of the system with colored 
multiplicative noise. Considering the noise amplitude as the power- 
law dependence of the stochastic variable x a we show the way to phase 
transitions disorder-order and order-disorder. The governed equa- 
tions for the order parameter and one-time correlator are obtained 
and investigated in details. The long-time asymptotes in the disor- 
dered and ordered domains on the phase portrait of the system are 
defined. 
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I.INTRODUCTION 

Considering an evolution of the stochastic system much more attention is fo- 
cused on the statistical moments behaviour. Usually, a system under consid- 
eration is described in terms of order parameter rj, being the first statistical 
moment (x), where x is a stochastic variable, and variance ((5x) 2 ), which 
plays role of autocorrelator (8x = x — (x)). Recently it has been shown 
that a multiplicative noise can lead a dynamical system to undergo a phase 
transition towards an ordered state (rj 7^ 0) fl], @. These results are obtained 
for extended systems within Curie-Weiss mean field theory and has been 
confirmed through extensive numerical simulations |§. In these works it was 
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shown that an ordered state exists only inside a window of noise intensities. 
Moreover, a spatial distributed system was investigated there and that kind 
of phase transition was stipulated by the symmetry breaking. The existence 
of the ordered state in the window of the noise intensity was explained by 
the collaboration between the multiplicative character of the noise and the 
presence of spatial coupling. Here we would like to describe peculiarities of 
the noise induced order-disorder phase transition for the zero-dimensional 
stochastic system. 

If we want to describe the evolution of our system we need to consider 
temporal dependencies of the averages according to average motion equation. 
More simplest case of ordinary thermodynamic system in this kind approach 
was described in [jjJJ, where a kinetic of phase transition is presented. A 
system with white multiplicative noise was discussed in ||, where we describe 
the evolution of the system with an arbitrary noise amplitude. Because the 
white noise is an idealized model of fluctuations of the system parameters, 
here we explore in what a way the colored noise can govern the system 
evolution. In this letter, instead of standard approach (small noise spectral 
parameter extension), we will use the unified colored noise approximation 
developed in || [7], |2|] to show the picture of phase transitions in stochastic 
systems. 

II.MODEL AND BASIC EQUATIONS 

In the simplest form, the problem of colored noise can be introduced con- 
sidering a relevant macrovariable x(t) (density of a given physical quantity) 
that satisfies a stochastic differential equation of the form 

x = f (x) + g (x)X(t), (1) 

where fo(x) represents a deterministic force. A stochastic part of the evolu- 
tion is defined by the amplitude g(x) of fluctuations introduced through the 
random term A(t), quite often assumed to be gaussianly distributed. With- 
out loss of generality, the deterministic part of the evolution can be chosen 
in polynomial Landau form 

/o(*) = -^> V{x) = - e -J + \x\ (2) 

where e is a parameter that acts as the dimensionless temperature, counted 
form the critical value. 
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Considering a whole set of models with a typical behaviour we can use 
the power-law function for the noise amplitude 



g (x)=x a , oe[0,l]. (3) 

This kind of assumption allows us to describe systems with self-similar phase 
space and in particular cases we can pass to the ordinary thermodynamic 
system (a = 0), directed percolation model (a = 1/2), population dynamics 
and forest fires (a = 1). 

In the simplest case for A(t) we can use definition of Ornstein-Uhlenbeck 
process 

rA = -A + £(t), (4) 
where r is the correlation time, is the white noise source ((£(£)) = 0, 

(£(*)«*)> = *(*-*))■ 

If we take the time derivative of Eq. ([!]), replace first A in terms of A and 
£ from Eq. (|4|) and then A in terms of x and x from Eq. (|l|) we can obtain the 
non-Markovian stochastic differential equation 

r(x-^x 2 ^ +ax = f + g o at), (5) 

where 

(6) 



a 
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According to the unified colored noise approximation we use the adiabatic 
elimination (neglecting x) and neglect x 2 . 

The problem now lies in obtaining an evolution equation for the order 
parameter and autocorrelator. For this purpose we average reduced equation 
of motion. It takes the form 

(a(x)x) = (f (x)). (7) 

The term in right hand side can be represented as a full derivative of dy(x) = 
adx, and after averaging, following ||, we get (dy(x)/dt) = (dy(x))/dt = 
(adx)/dt. Introducing the notation for the autocorrelator S = ((5x) 2 ) we 
rewrite right hand side of Eq. (0) as 

= T)(e + K(7] 2 + S)) + KT]S, (8) 
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where 



e = 1 — £t(1 — a), k = r(3 — a). 
The resulting equation for the first statistical moment reads 
[e + n(r] 2 + S))f] + nrjS = r](e — rj 2 ) — 3r]S. 



(9) 



(10) 



Because of Eq.(|H]) accounts effect of correlations we need to construct an 
equation for the autocorrelator. For this purpose let us express following 
differential: dy 2 = 2ydy + (dy) 2 . According to the aforementioned stochastic 
process dy(x(t)) we find dy 2 = e 2 dx 2 + (2en/3)dx 4 + (n/3) 2 dx 6 . Making use of 
the supposition i 6 C 1 we receive (dy 2 /dt) = e 2 d(x 2 )/dt + (2en/3)d(x 4 ) jdt. 
Rewriting the Langevin equation in the form of stochastic differential equa- 
tion 



dy = f (x)dt + g {x)dW 



we immediately produce up to the first order time derivative (where for Winer 
process we have (dW) 2 ~ dt) 



2er] 



e + 4k [ -rf + S 



ee{rf + S)-(e 



f] + e[e + 4k(t] 2 + S)]S = 
) (r] A + QrfS + 35 2 )1 + (x 2a ) 



(12) 



The obtained equation combines integer order averages and fractional 
one, namely (g^x)) = (x 2a ). As discussed in 0, such a kind of fractional 
average can be expressed in terms of order parameter and autocorrelator, 
according to the supposition that the probability distribution function of the 
initial process is a homogeneous function, i.e. 



P(x) w Ax' 



-2a 



A = -\1 - 2a\b ll - 2al . 
2 i I 



(13) 



where the cut-off parameter b 
the fractional average M: 



0. One yields the following definition for 

x nq ) = a n {q)(x n ) Pn{q) (14) 



4 



where 

1 — 2a + nq 

pM = t^t- (15) 

a n (q) = A^^p- l {q){l - 2a + nf"^ 1 . 

According to § a keypoint of the system with the multiplicative noise (§) is 
that its behaviour is governed by the magnitude of the exponent a in Eq.([|). 
Atl/2<a<l, when the fractal dimension of the phase space D = 2(1— a) is 
less than 1, the system is always disordered and its evolution is represented by 
the the autocorrelator S(t). It provides q > 1 in Eq . (|i4|) , hence the fractional 
order average is represented by S(t). In the case a < 1/2, where D > 1, 
according to the Landau theory the system can test the phase transition and 
in the fractional order average we account more essential contribution given 
by the order parameter rj (g < 1 in Eq.(|l4|)). Therefore, we can rewrite 
Eq.fll|) as 

,2a\ _ / "l^f 1 : < a < | 



where 



<^>H^' : i<a<l < 16 > 



2(1 ~ a) (17) 

III.EVOLUTION OF DISORDERED SYSTEM 

Let us analyze the evolution of the disordered system, at first. Considering 
the case a > 1/2 we get rj(t) = 0, hence the system evolution is defined by 
solutions of the following equation 

s(^ + 2kS) =S'(e-S , (3-— )) +a 2 S P2 . (18) 

The form of time dependencies for the autocorellator is shown in Fig.l. It is 
seen that S(t) monotonically attains the stationary magnitude determined 
by the equation 

— ) S + a^ 1 = (19) 



e 
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at condition Sq 7^ — e/4/t. In Fig. 2 we plot steady states at different values 
of noise correlation times and at different values of the control parameter 
e. According to Eq. (|i~9|) with e or r increase the stationary value So rises 
from the minimal magnitude. Let us focus on the limit S C 1. If we put 
5 P2 ^> 5 ^> 5 2 Eq. (|l8D gives the power-law time dependence 



t-»0 



Bt~^ B 



2(1 -p 2 )a 2 y 



(20) 



In the opposite case 5 — 5 <C 5 one has exponential dependence 5— 5 oc e xt , 
where A = (e — 5*0(3 — e/c/e) + a 2 P2)/(e/2 + 2nS ). 

IV.EVOLUTION OF ORDERING SYSTEM 



We now in a position to discuss the case a < 1/2 of ordering system with the 
fractal dimension in x — t space D > 1. It provides that the system dynamics 
is governed by equations for the order parameter and autocorrelator 



rj[e - rj 2 - 35] [e + 4/<7/ 2 + 5)] 



2K7] 
K€a\T\ 



e(rf + 5) 



KE 
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(r? 4 + 6r] 2 S + 35 2 



(21) 



P(v,s)s 



where 



= ~e( V 2 + S) - (l - + 6 V 2 S + 35 2 ) 



rj 



< l/v ( y + 5 



- ?7 

+ [e + «(77 2 + S)] ai rf 1 



[e - r? 2 - 35] 



7(^,5) = [e + ^ 2 + 5)][e + 4^ 2 + 5)] 

,2 



-2?fK 



r+ I/; ( |- + 5 



(22) 
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(3(r),S) = - + 2K(rf + S)\ [e + ntf + S)] 



(23) 



-KT) 



<■ ' i - ( y + S 



The obtained closed-loop system of differential equations can be ana- 
lyzed with a help of the phase plane method. From the corresponding phase 
portrait shown in Fig. 3a it is seen that at small values of the control pa- 
rameter e there is only one attractive point Co with coordinates i]q = 0, 
So = (e/3)[(l — £t(1 — o))/(l — 2er(l — 2a/3))]. Increasing of the control 
parameter provides the appearance of saddle and attractive points with co- 
ordinates given by solutions of the stationary equations 



[e - rj 2 -3S }[e + AK(rj 2 + S )} 



= 2k 



efaS + So) - (l - f ) (% 4 + H 2 S + 3S 2 ) 



(24) 



etfo + So) - (l - f ) (% 4 + 6% 2 S + 3S 2 )] x 
x [e + K{r,l + S )} = 

Sl [e- % 2 -3S ] 



(25) 



= Vl 



2 

' i j + So 



[e + k(^o + S'o)]ai77^ 1 . 



The dependence of the bifurcation value of the control parameter e from the 
noise correlation time r is shown in Fig. 4. It illustrates the appearance of the 
ordered phase domain (rjo ^ 0) at small values of the noise correlation time 
and small values of Eq. So the ordered phase is realized in the window of the 
intensities [£oj £ o]- Increasing of r provides the decreasing of the domain of 
ordered phase. Dependencies of the steady states are shown in Fig.5,6. Here 
thin lines display the saddle point S and thick lines correspond to the attrac- 
tive point C in Fig.3b. Some distinctive feature can be seen from Fig.5,6: 
the system undergoes the phase transition of the first order despite the bare 
a; 4 -potential corresponds to the continuous one. Such a kind of transition 
can be observed at disorder-order and order-disorder phase transitions with 
increasing and decreasing of the control parameter e. 

Let us discuss now time dependencies corresponding to the phase trajec- 
tories from Fig. 3. Time dependencies are shown in Fig. 7. 
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Thus, according to the phase portraits shown in Fig.5, at e ^ [^0' £ o] 
the order parameter falls down monotonically to the point Co, whereas the 
autocorrelator can vary nonmonotonically. Inside the domain [ejjj, £q], where 
bifurcation occurs, we can see two domains on the phase plane. These do- 
mains correspond to the small and large values of the order parameter. At 
small initial values of the order parameter we receive above mentioned be- 
haviour. At intermediate and large magnitudes the order parameter attains 
the attractive point C. In the vicinity of the saddle point S we get a critical 
slowing-down and metastable phase can exists some short time. Such a kind 
of behaviour can be found near the separatrix CqSC. 

Let us analyze the time dependencies of main averages analytically. Firstly, 
we investigate the system evolution in the disordered state (i.e. vicinity of 
the point Co) for large time. According to the || instead of the ordinary 
Lapunov method we use generalized Tsallis exponent JIU]: 



^ exp q (t) ^ [1 + (1 - q)^ (26) 

with generalized Lyapunov index q. This construction is more convenient 
to analyze the power-law dependencies that we have in evolution equations. 
According to the derivation rule 

^ exp,(t) = {exp q (t)) q = exp|(f) (27) 

and asymptotic behaviour 

limexp.(t) 1 + 1, lim expJt) -> [(1 - q)t] 1/1-9 (28) 



t—*oo 



let us assume solutions of Eqs. (|2l] , [22|) in the form 

rj(t) = mexp M (t), S(t) = S + nexp v (t). (29) 

Inserting Eq.(^) into Eq.(^T|) we receive up to the first order of m,n 1: 

(e + KS )(e + 4kS ) = -Kea im pi exp£ 1+1 ^(t), (30) 

where we account the singular contribution only. In the long-time limit, the 
function exp^ 1+1 ~ M (t) can be taken equal to 1. It provides definition of the 
Lyapunov multiplier and exponent 

m = \{e + KS ){e + AKS )/ne ai \ 1/pi , (31) 
1 



2(1 -a) 
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Considering behaviour of the autocorrelator S(t) up to the first order of 
amplitudes m and n we obtain 



(e/2 + 2kS ) = exp 1 -" (e + o^-W 1 exp^ 1 (t) exp; 1 (t)) . (32) 

In the short-time limit we assume exp* _t/ — > 1. The long-time asymptote 
yields exp^(t) exp~ 1 (t) = const = p^ 1 and provides 

zz = 2, n = — — . (33) 

nepi e/2 — e + ZKb 

Thus, according to the obtained time dependencies, we see that the order 
parameter behaves itself in a power-law form, i.e. r)(t) oc t~ 2 ^-~ a \ and the 
autocorrelator S(t) follows to the hyperbolical dependence S(t) oc t" 1 . 

If we pass through the critical value el then the system can be ordered and 
here we have to account an initial magnitude of the order parameter. Taking 
the initial magnitude of 77(0) larger then a critical value n c (shown in Fig. 8) 
we make the system to be ordered. Let us examine the system evolution 
to the ordered state (vicinity of the point C). Here we have to note, that 
we may not use the solution like generalized exponent (|26|). The latter is 
applicable at non-linearity effects which are sufficient to fix the amplitudes 
m and n. In the case under consideration the linear conditions are satisfied 
and because of power-law construction of the evolution equations we ought 
to use the Mellin transformations 



T)(t) = i] + J m q t q dq, 
S(t) = S + J n q t q dq. 



(34) 



Here the evolution equations can be transformed to the system of linear 
algebraical equations 

A n m q + A 12 n q = 0, 

a , a n ^35) 

A 21 m q + A 22 n q = 0, 

where multipliers Ay are functions of the system parameters and coordinates 
of the point C. The diagonal elements of the matrix A incorporate terms 
q/t, so we can rewrite An = A$ + q/t. The system ( p5| ) has solution when 



9 



det | A | = 0. So if we use the following notation c = —q/t we get 



T](t) = i] + mexp(— ct ln(t)), 
S(t) = 5*0 + nexp(— c£ln(t)), 



(36) 



here amplitudes m, n correspond to the index q = —ct; c is the real number 
whose magnitude can be expressed from the condition det |A| = 0: 



We have studied dynamics of noise induced order-disorder phase transition. 
Colored multiplicative noise contributes to realizing of two phase transitions 
with the control parameter increasing. We can conclude that this kind transi- 
tions are a consequence of collaboration of the nonlinearity of the system and 
multiplicative character of the colored noise (in the white noise limit there is 
only one phase transition along the axis e ||). We have shown that color of 
the noise does not change time asymptotes in the vicinity of the attractive 
points, corresponding to the disordered and ordered domains, comparing to 
the white noise limit. It changes only amplitudes of the time dependencies. 
We believe that these results may be found in physical and biological systems 
possessing a self-affine phase space. 
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FIGURE CAPTIONS 

Fig.l Time dependence of the autocorrelator S at e — 0.6, r = 0.5, 
a = 0.8 

Fig. 2 Stationary states of the system at a = 0.8: (a) So vs. control 
parameter at different values of r; (b) So vs. noise correlation time 
at different values of e. 

Fig.3 Phase portraits at a = 0.2: (a) e = 0.1, r = 0.01; (b) at e = 0.5, 
r = 0.2. 

Fig. 4 Phase diagram. 

Fig. 5 Stationary states of the system at a = 0.2: (a) order parameter 
vs. control parameter e; (b) order parameter vs. noise correlation 
time. 

Fig. 6 Stationary states of the system at a = 0.2: (a) autocorrelator vs. 
control parameter e (b) autocorrelator vs. noise correlation time. 

Fig. 7 Time dependencies correspond to the different trajectories on the 
phase portrait in Fig.3: (a) r\ vs. ln(i) at a = 0.2, e = 0.5, r = 0.2; 
(b) S vs. ln(t) at a = 0.2, e = 0.5, r = 0.2. 

Fig. 8 Critical value of the order parameter r\ c vs. noise correlation time 
t at e = 0.5 and different values of the noise exponent a. 
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